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1. INTRODUCTION 

Distance relay connected to a transmission line is required to operate for all types of fault within the 
protected zone. It is essential to determine correct type of fault in case of digital relaying for transmission line 
protection. The trip boundary is set according to the fault type as impedance seen by distance relay is not same 
for all types of fault [1]. Fault classification and location in a transmission line have been determined using 
modular artificial neural networks (ANNs) [2]. Optimal probabilistic neural network (PNN) with particle 
swarm optimization technique (PSO) has been considered for fault classification in a high voltage line [3]. 
PNN based method using Discrete wavelet transforms (DWT) are applied for classifying faults wherein DWT 
acts as an extractor of distinctive feature of input current signal [4]-[6]. Protection techniques based on 
different types of ANNs such as generalized neural networks, unsupervised/supervised neural networks, online 
trained neural networks and radial function basis neural networks are presented in [7]-[11]. Principle 
component analysis (PCA) based techniques have been applied to identify and classify faults in three phase 
transmission lines [12]. 

Literature shows neuro-fuzzy based techniques are widely considered for fault classification in 
modern power system involving complex networks [13]. ANN trained network can only be used in a system 
for which it is trained and implementing the same ANN to other network becomes cumbersome. An 8-bit 
microprocessor based technique for real time implementation is proposed in [14]. In this method voltage 
phasors are obtained using Kalman filtering and processed through a MC6809B microcontroller to distinguish 
between various classes of faults. The identification and classification of faults in HVDC system using wavelet 
analysis has been proposed in [15]. A deterministic fault classification method had been proposed in [16]. 
Linear perceptron based algorithm has been suggested to distinguish between grounded and ungrounded faults 
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for a distribution system [17]. The TMF has been used to differentiate between fault and a power swing to 
prevent maloperation of distance relay in zone-3 [18], [19]. 

This paper proposes a TMF based fault classifier wherein transient monitor (TM) indices of three 
phases are computed. The faulted phase TM index will exceed the threshold value and fault type can be 
identified through proposed logics. The phasor extraction of fundamental component of current signal is used 
to obtain the reconstructed signal. The phasor extraction accuracy of one-cycle DFT being high results in 
reduction of threshold value. The proposed method has been compared with wavelet transform based technique 
[16]. A 400 kV system has been simulated for different types of faults considering several factors. 


2. FAULT CLASSIFICATION TECHNIQUES 
2.1. Fault Classification Using DWT 

Discrete wavelet transform (DWT) can effectively analyse transient current and voltage. DWT 
analyse signals in different frequency band using short window for high frequency and long window for low 
frequency [20]. Therefore signal can be analysed in both time and frequency domain. Wavelet decomposes the 
signal in different frequency bands obtained from mother wavelet by dilation and translation. The wavelet 
transform of signal f(t) can be represented by following equation as: 


17 t-q 
WT(f, p.) =-= |f of) 
a P (1) 


where ‘p’, ‘q’ and ‘w’ are scaling constant, translation constant and wavelet function respectively. Discrete 
wavelet transform can be applied to transform of sampled waveforms as: 


1 


n-k : 
DWT(f,h,j)= grol 7 | 
Jp." P, (2) 


where ‘po!’ and ‘kpo” are scaling constant and translation constant respectively. 

The DWT can be implemented with successive pairs of high pass and low pass filters at each scaling 
stage. This is equivalent to successive approximation of same function where each approximation gives the 
incremental information referred to a particular scale. The shorter scale involves high frequency range and 
longer scale refers to low frequency of the spectrum. DWT decomposes a signal into a set of approximate and 
detailed coefficients. The wavelet transform based fault classification analysis involves mean of the 
approximate coefficients obtained from the DWT of each phase current at scale 1 and absolute value of the 
means are termed as Ma, Mp and M: respectively. Determination of fault type is carried out with logics shown 
in Table 1. 


Table 1. Fault Classification logics for DWT Technique 


Fault Type Fault Classification Logics 
A-G M, >My & M, > M. & M= M. 
B-G Mp >M, & Mp > M. & M= M. 
C-G M, >M, & M. > My & M= Mp 
A-B-G M, >M. & My > M, & Mat M> 2 M. 
B-C-G Mp >M, & M. > M, & Mit M. > 2 Ma 
C-A-G M, >My & M. > M, & M+M. 2> 2 My 
A-B M, >M, & My > M. & M= Me 
B-C Mp >M, & M. > M, & M= M. 
C-A Mp >M, & M, > Ms & M= M. 
A-B-C (Ma# Mp #M.)& (Ma + M, =M. )& (Mi + M. = Ma) & (Ma + M. = Mp) 


2.2. Transient Monitoring Function 

The fundamental component of current signal can be extracted by using a phasor estimation technique 
[1]. One cycle DFT has been considered for phasor estimation in order to reconstruct a signal from its phasor 
value at a given instant. In normal condition there will be no difference between the actual current signal and 
the reconstructed one. However any type of fault initiation results in significant difference between the two 
signals due to presence of transient component in fault signal. The absolute sum of squared difference between 
the actual signal and reconstructed signal over one cycle window length termed as TMF. Transient monitoring 
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had been applied to determine the accuracy of phasor estimation [20]. In this paper TMF is used as a fault 
classifier. The mathematical formulation of the TMF using Fourier transform approach has been outlined as 
follows. 

In the present paper fundamental frequency phasor value at each instant can be calculated using one 
cycle DFT. A data window of one cycle continuously moves forward and adds one sample corresponds to that 
instant. The phasor calculation for transient period will provide value for 50 Hz only. The current signal is 
reconstructed with phasor values obtained and compared with actual current signal. It can be observed from 
Figure 1 that during transient period two signals mismatches and TMF value increases. A sinusoidal current 
signal can be expressed as 


i, =I,, COS(ak+@) 


(3) 


where Im is the amplitude of current signal, k is the sample number, @ is the phase angle in radians and 
a=2m/N (N is one cycle window length). DFT of the sinusoidal current signal corresponding to fundamental 
frequency is 


Oodle 
l,=— Se oN 
NES (4) 
where Iı is the phasor estimated in terms of peak and can also be represented as 
= JP 
I, =1,,e (5) 
The signal can be reconstructed by using phasor values at each instant as 
i, = [L e Je"™ 
e A 


where T, is the sampling time and œo is the fundamental angular frequency. According to the definition of the 
TMF, it can be expressed as 


N 

TMF => C, 

= (7) 
where 

C, =i, —ix (8) 

At the inception of any type of fault, large transients will be observed in fault signal and TMF will 
maintain high index for few cycles. It can be observed from Figure 1 that actual current signal exactly matches 
with reconstructed current signal before fault initiation (at 0.5s). There occurs a mismatch during transient 
period although TMF maintains zero index value before fault and changing significantly during transient 
period. Fault classification is carried out using post fault current data of the three phases at one end of the line. 
Using these data TM indices for the three phases and ground are calculated and termed as TMa, TMp, TMce, and 
TMg where suffices a, b, c and g relate to phase A, phase B, phase C and ground respectively. 


TM, = ur SE + n) 


S 


(9) 
TM, -rmn Sv) 
Ts 
(10) 
TM, -rur 22 +x) 
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Where, FIT is the fault inception time. A delay of one cycle is considered for higher value of TM index. The 
total time taken for fault classification lies between one cycle and two cycles. The proposed fault classification 
technique requires accurate fault detection as a prerequisite. The logics for determining different types of faults 
are developed in terms of TM indices as illustrated in Table 2. 


Current,A 


Reconstructe 
i 1 


5 : A ; A 
(a) 0.46 0.47 0.48 0.49 05 0.51 052 0.53 


o 
(b) 0.46 0.47 048 049 05 0.51 0.52 0.53 
Time,s 


Figure 1. (a) Comparison between actual and reconstructed signal. (b) Transient monitor 


Table 2. Fault Classification logics for TMF method 
Fault Type Fault Classification Logics 


A-G TM, > h, & TM< h, & TM.< h, & TMe > hy 
B-G TM, < hp & TM,> hy & TM.< hy & TMs > hy 
C-G TM, < hy & TM< hp & TM. > hy & TM; > hg 
A-B-G TM, > hy & TM, > hy & TM< h, & TM; > hg 
B-C-G TM, < h, & TM,>h, & TM. > hy & TM, > he 
C-A-G TM, > hy & TMy< hy & TM. > hy & TM; > hg 
A-B TM, > h, & TMy> h, & TM. < hy & TM; < hg 
B-C TM, < h, & TMs > h, & TM. > hy & TM; < hg 
C-A TM, > h, & TMy< h, & TM. > hy & TM; < hg 
A-B-C TM, > hy & TM» > hy & TM. > hy & TMg < hg 


‘hp’ and ‘hg’ are the thresholds for phase and ground indexes respectively. 


3. RESULTS AND DISCUSSION 

A two terminal, 400 kV, 3 phases, 128 Km transmission line model has been considered as shown in 
Figure 2. The proposed classification method has been compared with DWT technique to show effectiveness 
and accuracy of TMF based classifier. All parameters related to the transmission line model are given in 
Appendix I. 


HV Transmission Line 


Fault 


Figure 2. Transmission Line model 


Extensive simulation studies using MATLAB®/SIMULINK™ have been carried out for variations in 
the values of fault resistance (Rp), fault inception angle (FIA), fault location (d) and load angle (ô). Normalized 
phase currents are processed and classification indices of DWT and TMF for phase and ground are obtained 
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for different types of fault under various system conditions. The indices value obtained for DWT and TMF 
method are given in Table 3 and Table 4 respectively. It is found from Table 3 that in few instances DWT 
indices fail to satisfy their corresponding logics with respect to fault type under particular system condition. 
Table 4 shows TMF indices satisfy their corresponding logics under the same system condition. 


Table 3. Values of DWT indices for various operating conditions 
Fault Condition: 


Fault Type d, Re, FIA, 8 [M] IMi] [M| IM,| 

0.5, 50Q, 20°, 20° 0.0024 0.0100 0.0069 0.0020 

A-G 0.1, 00Q, 45°, 20° 0.2353 0.0104 0.0273 0.0910 
0.9, 109, 25°,10° 0.0132 0.0121 0.0042 0.0018 

0.5, 50Q, 20°, 20° 0.002 4.334 6.691 0.571 

B-C-G 0.1, 05Q, 60°, 30° 0.071 76.73 121.1 3.848 
0.9, 100Q, 90°, 10° 0.006 1.470 1.600 0.113 

0.5, 10Q, 20°, 30° 0.000 13.53 13.24 0.000 

B-C 0.1, 02Q, 85°, 15° 0.000 727.3 731.8 0.000 
0.9, 50Q, 50°,20° 0.000 1.152 1.148 0.000 


Table 4. Values of TM indices for various operating conditions 
Fault Condition: 


Fault Type d, Ry, FIA, 8 TM, TM, TM. TM, 
0.5, 50Q, 20°, 20° 5.457 0.002 0.002 0.600 

A-G 0.1, 00Q, 45°, 20° 258.7 0.731 0.729 35.12 
0.9, 10Q, 25°,10° 49.68 0.267 0.267 4.028 

0.5, 50Q, 20°, 20° 0.002 4.334 6.691 0.571 

B-C-G 0.1, 05Q, 60°, 30° 0.071 76.73 121.1 3.848 
0.9, 100Q, 90°,10° 0.006 1.470 1.600 0.113 

0.5, 10Q, 20°, 30° 0.000 13.53 13.24 0.000 

B-C 0.1, 02Q, 85°, 15° 0.000 727.3 731.8 0.000 
0.9, 50Q, 50°,20° 0.000 1.152 1.148 0.000 


In literature different types of fault classification logics are found based on their fault classifier indices. 
The classifier index for three phases has been compared with each other. Alternatively the indices may be 
compared with a pre specified threshold value. In DWT technique |M| indices are compared with each other 
leading to inaccurate results due to a small margin between phase indices. It can be observed from Tables 1 
and 3 that logics based on comparison of phase indices results in false identification of fault type. TM indices 
are compared with a threshold value resulting in accurate identification of fault type. The value of thresholds 
‘hy’ and ‘hg’ are taken as 1 and 0.01 respectively. Optimal threshold is necessary to obtain accurate fault 
classification. Values of TM, and TM, for A-G fault are greater than hp and hg respectively whereas values of 
TM, and TM, smaller than hy. The results obtained shown in Table 4 validates the logics of Table 2 for other 
types of faults. 

Threshold value may depend upon constraints such as noise, harmonics, frequency variation and 
sampling rate. Phasor estimation quality affects the threshold value since phasor estimation accuracy decides 
similarity between actual and reconstructed signals. The difference between these two signals under prefault 
condition remains low leading to a very high value during fault conditions. One cycle delay after FIT is 
considered for high value of TM index as high threshold value avoids inaccurate classification. 

In the proposed method probability of correctness (Pe) is obtained for different cases to validate high 
classification accuracy. The parameter variation in each case over a long range has been studied in order to 
calculate P.. The sample space in each case consists of total number of parameter variation and favorable events 
pertaining to total number of correct classification. Mathematically ‘Pe’ can be expressed as: 


Number of correct classification 


Po = f f ; 100 (13) 
Number of iteration steps in each case 

The results using TMF shown in Table 5 reveal that the proposed method has high accuracy for wide 
range of operating conditions. The computation has been made for variation in Rr up to 130 Q, d and 6 up to 
0.95 (p.u) line length and 50° respectively. The transmission line resistance upto the fault point and fault 
resistance influence performance of proposed classifier. The low resistance results in high transient leading to 
increased mutual inductance between phases. High mutual inductance generate transient in healthy phases. TM 
indices of healthy phases increase due to presence of transients. The problem of high mutual inductance effect 
may be eliminated with the design of appropriate filters. It is found that proposed classifier is not affected due 
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variation in fault inception angle and load angle. The proposed technique is simple and straightforward as it 
requires a DFT based filtering algorithm and difference between two signals. 


Table 5. Classification accuracy of TMF method 


Fault Type Fixed Conditions Variable parameter Step Size Pe 
A-G d=0.5, FIA=45°, &=20° R= 00 to 1350 10 96.30 % 
B-C-G FIA=10°,6=10°,RF=50Q d=0.01 to 0.95 0.01 (p.u) 98.90 % 
A-C 5=10°,R=100Q, d=0.5 FIA = 0° to 360° 1° 100.00 % 

A-C-G Rr=1009,d=0.5,FIA=30° 5 = 5° to 50° 1° 100.0 


4. CONCLUSION 

A simple method for fault classification of power system has been proposed. The method employs a 
transient monitoring function based technique using current signals at one end of the line. Requisite logics with 
threshold values have been found to classify different type of fault. Simulation studies for various type of fault 
have been carried out. The effects of fault location, fault inception angle, fault resistance and load angle have 
been considered. It has been found that results obtained using proposed method are highly accurate compared 
to the Discrete Wavelet transform based method. 
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APPENDIX I 
The system information is given for power system operating at 50 Hz. For all instance negative 
sequence impedance is equal to positive sequence impedance. 
I. Source voltages: 
e Source 1: E;=400 kV 
e Source 2: Ex=40026 kV, where ô is the power angle. 
II. Source impedance: (Source 1) 
e Positive sequence impedance=1.744j19.92 Q 
e Zero sequence impedance=2.6100+j29.8858 Q 
IHI. Source impedance: (Source 2) 
e Positive sequence impedance=0.87+)9.96 Q 
e Zero sequence impedance=1.3050+j 14.9400 Q 
IV. Transmission line parameters (per km). 
e Positive sequence impedance=0.0256+j0.3670 Q 
e Zero sequence impedance=0.13624+j1.1096 Q 
V. Sampling Time: 1 millisecond. 
VI. Line Length: 128 km. 
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